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We discuss the one-dimensional, general quadratic Hamiltonian and the bi-dimensional charged 
particle in time-dependent electromagnetic fields through the Lie algebraic approach. Such method 
consists in finding a set of generators that form a closed Lie algebra in terms of which it is possible 
to express the Hamiltonian and the therefore the evolution operator. The evolution operator is then 
the starting point to obtain the propagator as well as the explicit form of the Heisenberg picture 
position and momentum operators. First, the set of generators forming a closed Lie algebra is 
identified for the general quadratic Hamiltonian. This algebra is later extended to study the the 
Hamiltonian of a charged particle in electromagnetic fields, given the similarities between the terms 
of these two Hamiltonians. 


I. INTRODUCTION 

The simple quantum oscillator is the building block of 
a very large number of well established physical models. 
Some of its most widespread applications are the atomic 
and molecular bonds that, under certain approximations, 
can be modelled by quadratic potentials. The time- 
dependent general harmonic oscillator (GHO), the most 
general version of a simple quantum harmonic oscillator, 
is at the heart of many interesting applications as radio¬ 
frequency ion traps. It consists of a simple harmonic 
oscillator with time-varying coefficients, time-dependent 
linear terms on the position and momentum operator and 
an extra term proportional to the symmetrized product 
of the position an momentum operators. It can be de¬ 
scribed by the quadratic Hamiltonian 

H = \a(t)p 2 + )^b (t) (xp + px) + ic (t) x 2 

+ d(t)p + e(t)x + g(t), (1) 

where x and p are the position and momentum opera¬ 
tors obeying the usual commutation relation \x,p\ = iTi 
and a, b , c, d, e and g are in general functions of time. 
Since in many cases it possesses exact solutions, it has 
turned into a key element to understanding and mod¬ 
elling a wide variety of physical systems where potentials 
are time-dependent. Specifically the GHO has been ap¬ 
plied in diverse branches of physics as quantum optics [1- 
3], transport theory in two dimensional electron systems 
[4, 5], quantum field theory [6], Ions traps (Paul traps) 
[7], laser cooling of trapped ions [8-10], quantum dissi¬ 
pation (Kanai-Caldirola Hamiltonians) [11-17], and even 
cosmology [18, 19]. One of the main advantages of mod¬ 
elling quantum physical systems with the GHO is that 
in many occasions it is exactly solvable [20]. The GHO 
has been studied by diverse mathematical methods such 
as the group-theoretical approach [21], the path integral 
approach [22], unitary transformations [7, 23], and the 


Lewis and Riesenfeld [24] invariant theory [16, 20, 25- 
28]. 

Besides the GHO, the time-dependent linear potential 
(LP), a particular case of the GHO, has also received con¬ 
siderable attention also due to the many applications in 
fields such as quantum optics, solid state physics, quan¬ 
tum field theory, molecular physics and quantum chem¬ 
istry among others. It has been established, at least since 
the 50’s, that the LP’s quantum propagator-and also the 
GHO’s propagator- possess a structure similar to the well 
known propagator for the simple quantum oscillator plus 
an interaction-dependent correction due to the forcing 
term in the Hamiltonian. [29, 30]. Whereas early studies 
of the LP relay on proposed Gaussian-like wave func¬ 
tion [31] and standard variables changes [32], recently, 
the quantum forced harmonic oscillator has been treated 
through more powerful methods as the Lewis and Riesen¬ 
feld [24] invariant theory [33-35], Feynman’s path inte¬ 
grals [29, 36-39], the generalization of the well known 
ladder operators [40], Laplace transform techniques [41] 
and time-space transformation methods [42]. 

Similarly to the GHO and the LP, the Hamiltonian 
describing a particle in time-dependent electromagnetic 
fields (CP) has countless applications in many physics 
fields such as quantum optics [43], single electron quan¬ 
tum dots [44] and magneto-transport theory [4, 5]. This 
system has been studied through different methods that 
include the Lewis and Riesenfeld [24] invariant theory 
[45, 46], path integral method [47], unitary transforma¬ 
tion approach [17, 48], and through quadratic invaritants 
[49]. 

The aim of this paper is to apply the Lie algebraic 
approach [50-54] to compute the evolution operator of 
the GHO. Drawing on these results we also calculate the 
evolution operator for the CP Hamiltonian. Additionally 
we obtain the propagator and the explicit form of the 
Heisenberg picture position and momentum operators. 

The mass-varying oscillator’s evolution operator was 
calculated by means of the SU (2) generators in Ref. [53, 


2 


55]. The Lie algebraic approach was also used to study 
the linear potential in Ref. [56] and the Kanai-Caldirola 
Hamiltonian in Ref. [57]. However, even though the Lie 
algebraic approach has been widely used to treat similar 
Hamiltonians, it has not been applied to solve specifically 
neither the CP nor the GHO Hamiltonians to the extent 
of our knowledge. 

The paper is organized as follows. In Sec. II we give 
an overview of the Lie algebraic approach. First, the 
time-dependent linear potential serves as an example to 
sketch the method and to work out some of the opera¬ 
tors that form the Lie algebra in Sec. III. Second, in 
Sec. IV, we deal with the evolution operator of the most 
general form of the quadratic Hamiltonian expanding the 
linear potential Lie algebra. With these general results 
we derive analytical expressions for a radio frequency ion 
trap in Sec. IV C and a Kanai-Caldirola forced harmonic 
oscillator in Sec. IV D. To complete our discussion we 
extend the Lie algebra of the GHO by introducing the 
angular momentum and extra generators. Finally we 
treat the Hamiltonian of a 2D charged particle in time- 
dependent electromagnetic fields exploiting the similari¬ 
ties of its Hamiltonian with GHO one. With the general 
expressions hereby obtained we compute analytical ex¬ 
pressions for a charged particle in time-varying magnetic 
field in Sec. V A and time-dependent electric fields in Sec. 
VB. We give the final conclusions in Sec. VI. The Lie al¬ 
gebra generators, their commutation relations and their 
structure constants are presented in Appendix A. Their 
corresponding unitary transformations are presented in 
Appendix B along with their transformation rules and 
propagators. 


II. OVERVIEW TO THE LIE ALGEBRAIC 
APPROACH 

A Hamiltonian is said to have a dynamical algebra if 
it can be expressed as the linear combination 

n 

H = J2ai{t)\u ( 2 ) 

2 = 1 

where ai ( t ) are real functions of time and the set of Her- 
mitian operators A = | Ai, A 2 ,... A n | forms a closed Lie 

algebra C. C is characterized by the structure constant 
c i,j,k hr the commutor 

n 

— ^ ^ ^ Afc. (3) 

fc =1 

In the sections to follow we show that the LP, the GHO 
and the CP Hamiltonians have dynamical algebras by 
identifying their generators and the corresponding struc¬ 
ture constants. 

The Lie algebraic approach [50-54] relays on the fact 
that the evolution operators of such Hamiltonians can be 


\i , A j 


expressed in either of the following forms 


U (t) = exp 


n 

i ^ ] OCi (t) A i , 
. 2=1 


U {t) = exp 

2=1 


iPi (t ) Aj 


( 4 ) 

( 5 ) 


where the transformation parameters at and @i are dif¬ 
ferentiable functions of time yet to be determined. 

We first consider Schrodinger’s equation 


H \il> (£)) = Pt | (t )), (6) 


where pt = ihd/dt , and, conveniently, we introduce the 
Floquet operator [58] 


H = H~Pu ( 7 ) 

that allows to write Schrodinger equation in the rather 

compact form 

H\1>(t))=0. ( 8 ) 

Using Eq. (5) let us now assume that there is a set of 
unitary transformations 

Q = {ui,u*, ■■■&„}, ( 9 ) 

with time-dependent transformation parameters /3i (t), 

/?2 (t), ■ ■ ■, P n ( t) in the form of Eq. (5) such that the 
application of 

u = u n ...u 2 Ui, ( 10 ) 

to the Floquet operator reduces it to the energy operator 
removing the Hamiltonian part as shown below 

uuW = -Pt. (11) 

We further assume that the explicit forms of the transfor¬ 
mation rules of Uixlfj , UipU} and UiptU.J for any unitary 
transformation in Q are known. The explicit form of the 
transformation rules of the unitary transformations used 
in this paper are presented in Appendix B. Conditions 
on the transformation parameters must be found so as 
to satisfy Eq. (11). As it is shown in Section IV, two 
different sets of unitary operators corresponding to the 
same Hamiltonian might comply with Eq. (11) meaning 
that there may be two or more different ways of arriving 
to the same evolution operator. 

If such a transformation does exist, the Schrodinger 
equation takes the form 

unWu (t)> = -p t (u |v> (t))) = 0 . ( 12 ) 

Reminding that pt is Ti times a time derivative, it is easy 
to see that U\tp (t)) is a constant ket, i. e. 

U (t)} = \tp (0)). 


(13) 
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Therefore the evolution of a quantum state ip can be 
easily calculated by multiplying the previous equation 
by the inverse of U (U _1 = W) getting 

\ip(t))=&\ipm. (i4) 

This equation states that obtained unitary operator U' 
is in fact the time evolution operator i. e. U* = U. 

The Green function, or the propagator, is calculated 
as usual in terms of the evolution operator as 


G (x, t; x', 0) = (x 


t/t 


X ) = { X 


= I dx 1 J dx 2 ... / dx n -l(x 
X (x 1 


X2 


ulul-.ul 

ul i 

Ul 


%n— 1 


X 

Xi 

x'). (15) 


The explicit form of the position and momentum opera¬ 
tors in the Heisenberg picture may be worked out from 
the transformation rules as 


x H ( : t ) = UxU\ 
Ph (t) = Upf7 t . 


(16) 

(17) 


A 4 = p 2 , where 1 is the identity operator. In appendix 
A 1 we present the commutors and structure constants 
of these generators; it is shown that in fact the algebra, 
exhibited by A 1 -A 4 , is closed. 

Even though this set of operators in principle guar¬ 
antees that the evolution operator should be given by 
U = exp (fa) exp (fax) exp (fap) exp (fa p 2 ), we proceed 
applying each generator’s unitary transformation step- 
wisely. Our first goal is to eliminate the linear term 
on the position operator x therefore we first apply the 
translation in space and momentum (see Appendix B 1) 
generated by Ai = 1, A 2 = x and A 3 = p given by 


Ui (t) = exp 


exp 

(t)x 

exp 

Y^(t)p 








where S (t), n (t) and A (t) are the real and differentiable 
time-dependent transformation parameters fa, fa and 
fa. The transformation rules for (20) are given by 

UiPtLl} = p t + S — An + tlx + Ap, (21) 

faxlll = x + A, (22) 

(hptll = p - n, (23) 


III. LINEAR POTENTIAL 


where an overdot denotes a time derivative. Under this 
transformation, the Floquet operator is transformed into 


To illustrate the use of the Lie algebraic approach, we 
analyze the solution of the one dimensional Schrodinger 
equation of a particle with variable mass subject to a 
time-dependent linear potential [33-35, 56]. The time- 
dependent mass term allows us to study dissipation in 
Kanai-Caldirola-like Hamiltonians [11-13]. Such Hamil¬ 
tonian is given by 

H = (18) 

where the mass m (t) and the force / (t) depend arbitrar¬ 
ily on time. From now on we drop their time-dependence 
except in special cases. The difficulty in finding the evo¬ 
lution operator for the time-dependent potential becomes 
evident when one computes the commutor of the Hamil¬ 
tonian at two different moments in times ii and £2 


H(h),H(t 2 ) 


ihp 


' f(h) 

_m(ti) 


/ 0 i) ~ 

m(t 2 )_ 


(19) 


UfaH - p t )Ut = ^-(p-H) 2 - f (t) (x + X) 

— (p t + s — An + ihc + Xp'j 

-(W-An-Q. (24) 

In order to vanish the linear terms in x and p we must set 
the following conditions on the transformation parame¬ 
ters 

n • 

- + A = 0, (25) 

m 

/ + n = 0, (26) 

with initial conditions A (0) = n (0) = 0 in order for U\ to 
be equal to the identity operator at t = 0 i. e. U\ (0) = 1. 
Equally, to cancel the independent terms, we must set 


In general, the last commutor does not vanish therefore 
the Hamiltonian (18) does not allow the evolution oper¬ 


ator to be written in the simple form exp 


-i f* H ( t ) dt 


requiring a different approach. 

Now we turn our attention to the generators of the 
Hamiltonian (18). At first glance, the set of operators 
{i,p 2 } seems like the right choice for £, however, a closer 
look at the commutation relations reveals that in order to 
close the algebra we must also include 1 and p. Thereby 
the whole set is given by Ai = 1, A 2 = x, A 3 = p, and 


n 2 

s + f A-An- —= 0, 

2m 


(27) 


with initial condition S (0) = 0. Immediately we notice 
the parallel between Eqs. (25) and (26) and the Hamil¬ 
ton equations of motion for the classical analog of (18). 
Moreover, if we collect the independent terms in Eq. (24) 
and define the classical Lagrangian 


^JL + An./A, 


(28) 
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its Euler equations yield the conditions imposed on the 
transformation parameters (25) and (26) 


Solving the differential equations (25)-(27) and (38) one 
obtains the transformation parameters 


d dL 
dt d\ 
d_dL 
dtdtl 


dL 

~d\ 

(LL 

dn 


n + / = o, 

(29) 

n 

-A = 0. 

m 

(30) 


The analogy goes even further when we notice that Eq. 
(27) is in fact the standard definition of the classical ac¬ 
tion 


A (i) = 

/ m (s) J dr/(r) ’ 

(41) 

n (t) = 

- [ dsf(s ), 

(42) 

Pit) = 

/'* ds 

Jo m (s) ’ 

(43) 


S = [ dsL(s). (31) 

J o 

Once these conditions are set, the original Floquet oper¬ 
ator is simplified into the one of a free particle 


Using Eqs. (16) and (17) we compute the position and 
momentum operator in the Heisenberg picture 

x H =x + fi(t)p+\(t), (44) 

Ph = P ~ n (f). (45) 


Ui {H - p t )U\ = J-p 2 -p t . (32) 

As it is desirable that all the terms from the Hamiltonian 
are eliminated, it is clear that the last transformation 
should be the one generated by A 4 = p 2 (see Appendix 
B 4) 


U 2 ( t ) = exp 


2h 


P (t) f 


(33) 


that yields the transformation rules 

Umul =p t + ^P 2 , 

U2XUI = X + P p, 

U 2 P &1 =p. 


(34) 

(35) 

(36) 


Application of this transformation to the Floquet opera¬ 
tor gives 


Hence, for given force / (t) and mass m ( t ) functions, one 
can easily determine all the transformation parameters 
through Eqs. (41)-(43) and plug this solutions into the 
propagator and Heisenberg picture space and momentum 
operators. 

Finally, from Eq. (15) and the propagators shown in 
Appendices B 1 and B 4, the propagator for the LP can 
be expressed as 


G (x, f; x', 0) = / dx\ (x Ul 


x\ ) ( X\ 


m 


T%P ( t ) 


exp 




n (t) 

exp ^ — 1 —-— [x — A ( t)\ 


exp ^ [x — x' — A (t )] 2 . (46) 


2hp ( t ) 


IV. GENERAL QUADRATIC HAMILTONIAN 


U 2 U 1 (H - pt) U\UI =\[y i ~ti)p 2 - Pt- (37) 

By establishing the restriction 

--4 = 0 , (38) 

m 

with the initial condition fd ( 0 ) =0 in order to make 
U 2 (0) = 1, the Floquet operator is finally reduced to the 
energy operator 

U 2 U^H-pPiUlu^-pt. (39) 


With the earlier treatment we can handle the GHO 
Hamiltonian with time-dependent coefficients. We follow 
two different procedures to obtain the unitary operator 
for the GHO Hamiltonian in order to study two different 
special cases: a radio frequency ion trap and the Kanai- 
Caldirola Hamiltonian of a forced harmonic oscillator. 

We start with the most general Hamiltonian 

H = \a (t) p 2 + h) (t) (xp +px) + ^c ( t) x 2 

+ d(t)p + e{t)x + g(t), (47) 


Therefore, by Eqs. (13) and (14) the evolution operator 
is given by 


(£) = u\ (£) ul (£) = exp 




x exp 


~^A (t)p 


exp 

x exp 


-~n{t)x 




(40) 


with arbitrary time-dependent coefficients a(t), 6 (t), 
c(t), d(t), e(t) and g(t). It is quite clear from the 
structure of (47) that the closed algebra corresponding to 
this Hamiltonian should be given by the set of operators 
Ai = i, A 2 = x, A 3 = p, A 4 = x 2 , A 5 = p 2 , A 6 = xp + px. 
In Appendix A 2 we present the commutation relations 
and the structure constants for these generators. 

We first aim to remove the independent terms and the 
ones proportional to x and p. We thus apply the trans¬ 
lation in space and momentum shown in Appendix B 1 
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generated by Ai = 1, A 2 = x and A 3 = p given by 


U\ (t) = exp 


exp 

in(t)i 

exp 







( 


The transformation rules for (48) are given by 

UiPt U\ = Pt + S - An + tlx + A p, (49) 

frixUl = x + A, (50) 

UxpUl =p- n. (51) 

Applying (49)-(51) to (47) The transformed Floquet op¬ 
erator takes the form 


Ui (h - p^j u\ = i ap 2 + (xp + px) + ^ cx 2 - p t 
+ x ^cA — 611 + e — II) 

T p ^—all T bX T d — a) 

+ g - S + ian 2 + icA 2 - 6An - dU + eX + AIL (52) 

It is possible vanish the independent and linear terms 
in x and p by imposing the following restriction on the 
transformation parameters 

cA - bll + e - II = 0, (53) 

-an + bX + d - A = 0, (54) 


S = g + —all 2 + ^cA 2 — 6An — c?n + eA + An, (55) 

with initial conditions S (0) = A (0) = n (0) = 0 in order 
to guarantee that U\ (0) is the identity operator at t = 0. 

Once again, we observe that Eqs. (53)-(54) are 
the classical Euler equations corresponding to the La- 
grangian 

L = ian 2 + tjcA 2 — 6 An — cdl + eA + An + g. (56) 


We can readily obtain Eqs. (53) and (54) from the Euler 
equations 


d dL dL 

dt dX dX 

d dL dL 

dt dtl d\\ 


ri — cA + fen — e = 0, 
A + an — bX — d = 0, 


(57) 

(58) 


and Eq. (55) yields the very well known relation for the 
action 


S = 



(59) 


Imposing conditions (53)-(55), the transformed Flo¬ 
quet operator reduces to the quadratic form 


Ui (h - U\ 


= ]^ap 2 + ^b(xp + px) + 7 jCX 2 - p t . (60) 


As pointed out earlier, there may be different sets of 
unitary transformations that reduce the Floquet operator 
to pt- To illustrate two possible solutions for the evolu¬ 
tion operator of the general quadratic Hamiltonian, at 
this point, we take two different calculation paths. The 
first one is shorter but requires the solution of Riccati 
differential equation whereas the second path is more in¬ 
volved but in a wide variety of physical situations avoids 
solving Riccati differential equation through the use of 
Arnold transformation. 

Now we consider the dilation generated by Ag = xp+px 
(see Appendix B 2) given by 


U 2 (t) = exp 


2 k 


7 (t) (xp + px) 


that yields the following transformation rules 

U 2 p t U 2 =Pt + ^7 (xp + px ), 
U 2 xU\ = e 7 x, 

U 2 pU\ = e-Tp. 


(61) 


(62) 

(63) 

(64) 


The application of the dilation yields the following trans¬ 
formed Floquet operator 


U2U 1 (H-Pt) Ul&l = ^ae“ 2 V + \c^x 2 

+ \ 0- 7 ) (xp+px) -Pt- (65) 

Although doing 7 = b to remove the term proportional to 
(xp + px) would seem to simplify the Floquet operator, 
it is more convenient to leave 7 as a free parameter that 
will be useful later on. 


A. First path 


We take the calculation from Eq. (65). In this path it 
is convenient to set the 7 parameter by doing 

e 27 = aA, ( 66 ) 

and A = 1/a (0) = 1/ao in order to make the dilation U 2 
equal to the identity operator at t = 0. This seemingly 
arbitrary definition of 7 and A will prove to be a key step 
in simplifying Riccati equation into a linear second order 
differential equation. 

It is possible to get a notable simplification by applying 
the unitary transformation generated by A 4 = x 2 (see 
Appendix B 3) given by 


U 3 = exp 


2 h 


a (t) Ax 2 


that yields the following transformation rules 

A 2 

U-iPtU^ =p t +a—x , 

U 3 xUl = x, 

U 3 pU\ = p — aAx. 


(67) 


( 68 ) 

(69) 

(70) 
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After carrying the transformations above the Floquet op- In order to eliminate the term proportional to xp + px 
erator takes the form we set 


U3&2U1 ( H-pt) U\u\ut 

= — [ae~ 27 A 2 a 2 — 2 (6 — 7 ) Aa + ce 27 — Ad] x 2 

+ -ae -27 p 2 + — (b — 7 — ae _ 27 Aa) ( xp + px) — pt- 

(71) 

It is desirable that the term proportional to x 2 vanish, 
hence we restrict the values of the a parameter by setting 
the condition 

Ad = ae~ 2y A 2 a 2 — 2 (b — 7 ) Aa + ce 27 . (72) 


cj) = b — 7 — ae 27 Aa, (81) 

obtaining the Floquet operator for a free particle with 
variable mass 

U4U 3 u 2 Ui (A-&) ii\ulu}M 

= i ae - 2 ( 7+ « P 2 - pt. (82) 

It is clear that in order to remove the remaining term 
we must apply the transformation generated by A 5 = p 2 
that can be expressed as follows 


This is a Riccati differential equation of the form 

y' (x) = Qo (x) + qi (x) y (z) + q 2 (x) y 2 (x) , (73) 

with a = y, qo = ce 27 /A, qi = — 2(6 — 7 ) an d <?2 = 
ae _ 27 A. Note that applying the restriction ( 66 ) and do¬ 
ing the variable change a = — uju , the Riccati equation 
is turned into the simpler linear second order differential 
equation 


u + 



u + cau = 0 . 


(74) 


If Eq. (72)- or equivalently Eq. (74)- hold the trans¬ 
formed Floquet operator becomes 


U 3 U 2 U 1 (H - p t ) Ululilt = \ae ~ 27 f 


+ — (b — 7 — ae 27 Aa) (xp + px) — pt- (75) 
The dilation generated by Xq = xp + px 


U 5 (t) = exp ^r/3 (t) ^ , 

(83) 

with the transformation rules given by 


UbPtUt =Pt + /3-^p 2 , 

(84) 

U 5 X&I =x + /3^, 

(85) 

U 5 P&I = P- 

( 86 ) 


Therefore, under this transformation the Floquet opera¬ 
tor takes the form 

U5U4U3U2U1 ( H-pt) uMululut 


1 

2 


ae -2( 7 +0) _ 


£ 

A 


Pt- (87) 


Imposing the following restriction on the transformation 
parameter 


Ua ( t) = exp 




-px) 


(76) 


seems the right choice for the next transformation since 
it trivially commutes with its generator xp + px. This 
transformation produces the transformation rules given 

by 


UapM =Pt + ^4> (xp + px ), 

(77) 

UaxUI = e^x, 

(78) 

Ua0\ = e_ 0 P. 

(79) 


Application of the dilation yields the transformed Flo¬ 
quet operator 


UAU3U2U1 ( H-pt) ulululul = iae- 2 ( 7+ « f 


+ i (b — 7 — <j) — ae 27 Aa) (xp + px) — p t . (80) 


13 = A ae- 2(7+ ^, ( 88 ) 

the Floquet operator is finally rendered into the energy 
operator i. e. 

U 5 U 4 .U 3 U 2 U 1 (H - p t ) Ultjluiului = -p t . (89) 

According to Eqs. (12) and (13) this product of trans¬ 
formations is precisely the evolution operator 


C/t = ululululul = exp ( — yAp ) exp ( -^ITr 


x exp ( ——S ) exp 
exp ( ——aAx 2 ) exp 


~2jA(xp + px) 
-^rfiixp+px) 


X6XP (90) 
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Collecting the results above, the classical equations of 
motion for the position and momentum are obtained from 
the Euler equations (53) and (54) 


Using Eq. (15) and the Green functions of the five uni¬ 
tary transformations in Appendix B we readily integrate 
the propagator 


li = cA - MI + e, (91) 

A = bX - an + d, (92) 

and the corresponding classical action can be calculated 
by substituting the explicit forms of A and II into (55) 
and integrating 


S = 



9 0) + \ a 0) n 2 0) + 0) A (s) 2 

- b(s) \{s)H(s) - d(s)U(s) 


+e ( s ) A (s) + A (s) II (s) 


(93) 


To complete the remaining parameters, from ( 66 ) we set 
7 = ln(a/ao)/2 in order to simplify the Riccati Eq. (72). 
Next, the a parameter may be integrated either from (72) 
or (74). Then, the <j> and f3 parameters are calculated by 
direct integration of the ordinary differential equations 
(81) and ( 88 ) 


G ( x , t\ x', 0) = 
x exp 


2irh/3 

Ae “ 27 fe~ 24. 


exp 


— i- 


: S{t ) 


exp - 


^ + 7 


2 n 

x exp 


P 


— a (x — A)"* 


exp 


■ A 2 

O ) 


2h/3 


“• {—hp~ x + a ) (I “ A) 


( 102 ) 


B. Second path 

Here we follow an alternative path to the one in the 
previous section. In this path we start the calculation 
from Eq. (65) but instead of ( 66 ) we impose the following 
restriction on the 7 parameter 


(103) 



Pit) 


-7 it) 


+ ds 


b(s) — a (s) e~ 27(s) Aa (s) 


dsAa (s) e-^W+Ml. 


(94) 

(95) 


Putting the explicit form of the evolution operator (90) 
into Eqs. (16), (17) the Heisenberg picture position and 
momentum operators can be expressed as 


where A = tJcq/oq in order to make the dilation U 2 
equal to the identity operator at t = 0 . 

First we apply Arnold’s transformation 


U 3 (t) = exp 




(104) 


where <f> is the transformation parameter. The transfor¬ 
mation rules for (104) are 


where 


and 


x H it) 

PH it) 



A 

-n ’ 


m = 


G qq G qp 

t~*pq Gpp 


a = e ^+ 7 

^ qq — c j 

r — 7 

'-’IP ~ ! 

G pq = -aAe*~\ 

Gpp = e-*-' 1 - ape^. 


(96) 


(97) 


(98) 

(99) 

( 100 ) 

( 101 ) 


The matrix M satisfies the symplectic conditions in¬ 
herited from the unitary transformations M T icr y M = 
MicTyM 1 = itjy with (T y the Pauli matrix. Additionally 
it complies with detM = 1. The three previous con¬ 
ditions ensure that the commutation relations between 
position and momentum operators are preserved during 
the system’s evolution, namely [x#(t),p# (t)] = ih. 


Umut =Pt + ^<j> + Ax 2 ^j, 

(105) 

U 3 xUl = x cos p + —p sin </>, 

(106) 

U 3 pUl = p cos p — Ax sin </>. 

(107) 


Note that Arnold’s transformation is generated by a lin¬ 
ear combination of A4 = x 2 and A5 = p 2 . Under this 
transformation the Floquet operator takes the form 


U 3 &2Ui (H-Pt) U\ulU: 


t _ 
3 — 


+ 


1 

2 A L 
A ■ 
2 . 


ac — <j> + (b — 7 ) sin 2<p 


ac — 4> — {b — 7 ) sin 2 <j) 


+ ^ (6 - 7 ) cos 20 (xp + pi) —pt- (108) 

By restricting Arnold’s transformation parameter by 
the relation 


1 = y/ac , 


(109) 






































the Floquet operator reduces to the following quadratic 
form proportional to b — 7 


8 

After the condition (118) has been set, the Floquet 
operator (117) takes the form 


U 3 U 2 Ui ( H-Pt) ututllr 




A 


— Ax 2 I sin 20 


Certain cases where 


+ (xp + px ) cos 20 


b- 7 = 0, 


-Vt- (no) 


( 111 ) 


U 4 U 3 U 2 lh ( H-Pt) ululutil] 


= 2 (b-V 


V P 2 


sin 20 + (xp + px) (cos 20 — a sin 20) 


-Pt- (H9) 


It is convenient to set the nearly last transformation 
to be a dilation of the form 


U 5 ( t ) = exp 


2 h 


ip (t) (xp + px) 


( 120 ) 


specially those where b = 7 = 0, lead to physically mean¬ 
ingful systems such as a variable mass charged particle 
in constant magnetic field. It is thus worthwhile to treat 
them separately. If condition (111) is fulfilled the Floquet 
operator is completely reduced to the energy operator pt 
implying that the evolution operator is simply given by 

u' = u\utui ( 112 ) 

However, in order to consider cases where b — 7 0 0 we 
must move on to the next transformation. We consider 
the unitary transformation generated by A4 = x 2 given 

by 


U 4 (t) = exp 


±a( t )A^ 


(113) 


since it just multiplies the p 2 term by a factor exp (—2<p) 
and yields an additional ip (xp + px) term that allows to 
cancel the (b — 7) (xp + px) (cos 20 — a sin 20). Indeed, 
transcribing the transformation rules from (62)-(64) 


Usptty =Pt+ ^ A {xp + px ), 

(121) 

U 5 xUl = e^x, 

(122) 

U 5 pUl = e“^p, 

(123) 


and applying t/5 to the previous Floquet operator we 
obtain the above described terms 

U 6 U 4 U 3 U 2 Ui (H-Pt) uimtulut 

= ^ (b - 7) e~ 2v p 2 sin 20 


with the following transformation rules 


M A 2 

UiPtUl =p t + a—x , 

(114) 

U 4 xU\ = x, 

(115) 

U 4 pU\ = p — a Ax. 

(116) 


The application of this transformation to the Floquet 
operator in Eq. (110) yields 


+ 2 + P&) [(& — 7) (cos 20 — a sin 20) — p] 

-Pt- (124) 

By imposing 

p = (b 7) (cos 20 gl sin 20), (125) 

the Floquet operator is reduced to the one of a free par¬ 
ticle with variable mass 


UiU 3 U 2 Ui (H-Pt) ulululu) 


= - 2 (b-i) 


— sin 20 + (xp + px) (cos 20 — a sin 20) 


+ - {(b — 7) [(a 2 — l) sin20 — 2cxcos20] — d} Ax 2 . 


(117) 


To eliminate the terms proportional to x 2 we set the fol¬ 
lowing restriction on a 

a = (b — 7) [(a 2 — l) sin 20 — 2a cos 20] . (118) 

This is newly a Riccati differential equation of the form 
(73) with q 0 = — (b — 7) sin 20, q 3 = 2 (b — 7) cos 20 and 
<72 = (b - 7) sin 20 . 


U 5 U 4 U 3 U 2 Ui (H-Pt) utmtulut 

= ^ (6-7)e“ 2v p 2 si n2 0-j5 t . (126) 

Evidently, the last transformation to be used is the one 
generated by A5 = p 2 


U 6 (t) = exp ^-/3 (t) ^ , 

(127) 

with transformation rules 

U&PtUl = p t + P-^p 2 , 

(128) 

UqxUq =x + P^, 

(129) 

U 6 P&e = P- 

(130) 
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In this case, the Floquet operator takes the form 


uMMsUi ( H-pt ) ulu^ulu^ 


1 r 


= (6 — 7 ) e 2<p sin 20 /3 -- - p t - (131) 


2 L 


The Floquet operator above can be reduced to the energy 
operator p t by imposing the following condition on /3 

/?=(& — 7 ) e~ 2lp sin 2(j). (132) 

We have thus arrived to the form (11) and therefore by 
collecting (48), (61), (104), (113), (120) and (127) the 
evolution operator is given by 


m = ululutuluiul 

= exp ( — — S ) exp ( — — Xp ) exp ( ——ni 


h 

x exp 


-^7 (xp + px) 


x exp 


( A ^ 2 + TP 


x exp 


2 h 


(xp + px) 


exp ( ——aAx 2 


exp 


■—/3— 

2 rA 


(133) 


The transformation parameters must be calculated from 
the restrictions (53), (54), (55), (103), (109), (118), (125) 
and (132). 

By successively applying the six transformations above 
to x and p we can workout the Heisenberg picture posi¬ 
tion and momentum operators newly obtaining the sym- 
plectic form 


where 


x H (t) 

PH (t) 



A 

-n ’ 


m = 


Gqq Gqp 
Gpq Gpp 


(134) 


(135) 


but instead, in this case the matrix elements are given 

by 


Gqq (t) = (cos (j) — a sin </) e y+ip 


= (cos (j) — a sin <f>) (136) 

a 0 c 


Gqp (t) = [(/3 cos <j) — a(3 sin cj>) e v + sin </e ^ 

= I "(P cos <t> — ot/3 sin (j>) e v + sin <^e _¥, l , (137) 

V c 0 c 


Gpq (t) = — (a cos (j> + sin (j>) Ae^ 7 


= — (acos</> + sin</>) (138) 

OqU 


Gpp (t ) = — [(/3 sin (/) + a/3 cos </) e v — cos fie 

= — [(/3 sin (f> + a/3 cos (j>) e^ — cos (j)e~ v 1 , 4 /— 

V Coa 


(139) 


Even though the structure of the matrix M is radically 
different from the one obtained in the first path [see 
Eqs. (98)-(101)], it also satisfies the symplectic condi¬ 
tions. From Eq. (15) and the propagators presented in 
Appendix B we can obtain the propagator associated to 
(133) as 


G 


(x,t\ x',0) = J 


= / dxidx2dx3dx4dx5 


X ( X 




Xi) { X\ 


Ui 


X ... < x 5 


Ut 


X2j [X 2 
X 6 


A 2 


4iTrh 2 f3l sin <5 


.S (p + 7 
6XP 1 -'h - ~2~ 


n 


ui 


exp 


+iu (pc') + i ( qx' — — ) (a: — A) 


X 3 

iw (x — A ) 2 
, (140) 


where, for the sake of brevity, we have defined the fol¬ 
lowing functions 


A / Ae-^A 
“ “ Zhj3 V + 2 h/3l ) ’ 

Ae~ 2ip f A \ 

2 /ism (p \ 2hlsm(f)J 

A2 e -(<P+7) 

q ~ 2h 2 l(3 sin^’ 

, A /3 (a — cot d>) — e~ 2v> 

l = A— - — -. 

2/i/3 


(141) 

(142) 

(143) 

(144) 


C. Radio frequency ion trap 

Here we assume that the trap potential can be decom¬ 
posed into a static and a time-dependent part that varies 
sinusoidally at the drive radio-frequency to [9]. Thus, the 
ion Hamiltonian is given by 

H = + —A -(- -)- — (K x + k x cos wt) x 2 

2m 2m 2m 2 

+ — ( K y + k y cos lot) y 2 + — ( I\ z + k z cos ujt) z 2 . (145) 

The parameters k x , k yi k Zl K x , K y and K z are restricted 
being that the electrical potential has to fulfil Laplace 
equation. A possible choice is — ( K x + K y ) = K z > 0 
and k x = —k y . The separability of the previous Hamil¬ 
tonian allows us to work the x , y and z coordinates inde¬ 
pendently. Thereby cq = 1/m, b z = 0, Ci = Ki + kt cos wt, 
di = ei = gi = 0 where i = x,y,z. The solution to the 
ordinary differential equations (91) and (92) with initial 
conditions (0) = A j (0) is n, = A; = 0 and thus, the 
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action is S = 0. When 92 = 1 Riccati Eq. (73) is known 
to reduce to a second order linear equation by making 
at = —iii/ui . Therefore, we set Aa^ exp (—27,) = 1 by 
doing Aj = m and 7* = 0. We are left with the Mathieu 
differential equation 


Ui + — (Ki + ki cos oot) m = 0, (146) 


whose solution is given by 

4 Ki 


t, = AC 


mui* 


2 ki Lot 

2’T 


mui 


(147) 


where A is a constant and C (a, q , z ) is the Mathieu co¬ 
sine function that complies with C (a, q, 0) = 1 and its 
derivative C (a,q, 0) = 0. The parameter is therefore 
given by 


V rriuj ^mar 2 / 

\ muj^ 1 mui* 1 1 / 

By substituting the previous result in (94) we get 


(148) 


<j>i = - [ dsoti ( s ) = [ 

J o Jo Ui (s) 


= In 


C 


4 Ki 


2k j cut 


rriLd z 


■ (149) 


The last parameter is 


D. Forced harmonic oscillator with varying mass 


The harmonic oscillator with varying mass Hamilto¬ 
nian is a useful theoretical tool to study quantum dis¬ 
sipation [11, 12]. It has been treated by diverse meth¬ 
ods including Feynman integrals [23], and the Lie alge¬ 
braic approach [53]. In the latter the Hamiltonian was 
expressed by means of the three SU(2) generators. In 
contrast, the forced harmonic oscillator needs a larger 
set of generators because of the linear potential terms. 
The expression for the Kanai-Caldirola Hamiltonian of 
the forced harmonic oscillator is 


H = 


=-t/r 
2m 1 


ct/r 


2 ~2 


-mu n : v 


— e 4 / T (Fq + F\ sin wit) x. (154) 


Therefore, a = exp {—t/r) /to, c = Wgmexp (t/r), e = 
— exp (t/r) (Fq + F\ sinuqt) and b = d = g = 0. Re¬ 
stricting ourselves to the case of over-damping, i. e. 
4 t 2 Wq < 1, the A and n parameters can be obtained 
from the solution of the ordinary differential Eqs. (91) 
and (92) that yield the standard solutions for the classical 
damped harmonic oscillator 


A(f) 


Fq 

muiQ 


^1 — e */ 2r coshnt 


e-*/2r 
2 rfi 


sinh f It 




Pi = T 


4 Ki 


2 ki cot 
moo 2 ’ 2 


\ /•' ds 

4- 

Fi 

) I / 4 Ki 2 kj us\ 

^ ® V muj' 2 ’ muj' 2 ’ 2 / 

m 

T 2 (uf-Lotf+CO 2 


(150) 

Gathering the results above, the Heisenberg picture op¬ 
erators are given by 

^(t) = c i (f)f + ic i (f)r i (f)A,( 151) 


PiH it) = 




Pi 


moo / Lot\ „ 

+ ~ Ci \Y ) Xh 


(152) 


TLO\e t d 2T cosh f It + —(l + 2 t 2 uj\ — 2 t 2 Wq) e t ^ T 
x sinh Qt — tloi cos u>\t + t 2 (w 2 — oof) sin oj-pt , (155) 


H{t) = —jfe t/2T sinh nt + 


T-i 


r 2 (lo 2 ~lo 2 ) 2 +u; 2 


t 2 lo± (wq — w 2 ) e 4 / 2r ^ cosh AM — e 4 / 2r cos uoitj 


where, for the sake of simplicity we have de¬ 
fined Ci(ujt/ 2) = C (AKi/muo 2 , — 2k t / muo 2 , uot / 2) and 

Ti (oot/2) = fg dsC~ 2 (AKi/mu) 2 ,—2ki/mui 2 ,ws/2). Fi¬ 
nally, the propagator is given by 


G{x,y,z,t\x',y',z', 0) = 



et/,2T sinh fit — Tu(e 4 ^ sinwif 


(156) 


where H = Y\ \ — 4 t 2 Wq|/ 2 t. 

Now we turn to the a, <j> and /? parameters. As in the 
previous example, (72) is the key equation we have to 
solve first. To do so, we set 7 = —t/2r, A = m and do 
the variable change a = —u/u rendering Riccati equation 
in the form of a damped harmonic oscillator 

1 9 , 

Li -I —u + loqU = 0. (157) 

T 


Under over-damping conditions, the solution for the pre¬ 
vious differential equation is 


a(t) 


1 1 - 4r 2 H 2 

2r 1 + 2 tQ, coth f It 


(158) 
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Integrating (94) and (95) we obtain the remaining pa¬ 
rameters 


cf) (t) = In 


^cosh fit + 


1 

2rfi 


sinh fit 


m 


2 T 

1 + 2 rfl coth fit 


(159) 

(160) 


The under-damped harmonic oscillator parameters are 
obtained by doing fl —> ifl in (155), (156), (158), (159) 
and (160). Note that the three previous results are com¬ 
parable to the ones obtained in Ref. [53] by using the 
SU( 2) generators. Substituting the explicit forms of the 
parameters into (96) we obtain the Heisenberg picture 
position and momentum operators 


where x, y, p x and p y are the standard space and mo¬ 
mentum operators in the x — y plane. The electron’s 
charge is given by e and the scalar and vector potentials 
are expressed in the completely symmetric gauge by 


E x ( t) x- Ey (t) y, 

(165) 


(166) 

B (t) x. 

(167) 


The mass to, the magnetic field B and the coefficient 
K may be time-dependent. Substituting the scalar and 
vector potentials in the expression for the Hamiltonian 
and expanding, we obtain [45] 


x H (t) 


( cosh fit + —— sinh fit ) e t ^ 2j 


2 rfl 


e -t/2T 

H-—- sinh fit p + A, (161) 

mil 

Ph (t) = (4r 2 fl 2 — l) e i / 2r sinh fit x 

4 T Z ll 

e t/2r 

-\ -— (2 rfl cosh fit — sinh fit) p — n.(162) 

2 rfl 

Finally, introducing the explicit form of the parameters 
into (102) the propagator can be expressed as 


G (x, t; x', 0) = 


mfl 


2nh sinh fit 


e ~iS/h gt/4r 


x exp 

x exp 


—*—— e^ r (1 — 2rfl coth fit) (x — XY 
Ahr 


TYl o' 

i —— (1 + 2 rfl coth fit) ( x ') 


. Ahr 


exp H(isSm e ‘ /2Tl ' + 3 <3:_A) ]' (163) 


The Heisenberg picture position and momentum oper¬ 
ators and the propagator in the under-damping regime 
can easily be found by doing the H —> ifl. 




eB 

2m 


(xp y 



- VPx) + eE x x + eE v y. 


(168) 


The structure shown by this Hamiltonian (168) suggests 
that the set of generators that yields the corresponding 
closed Lie algebra should be at least composed of the 
identity operator, the generators listed in the previous 
section (A2 to Ae) for the x and y parts of (168) and 
the angular momentum L z = xp y — yp x . However three 
more generators are needed in order to close the algebra: 
xp y + yp x , xy and p x p y . Thereby, the complete set is 
given by Ai = 1, A 2 = x, A 3 = p x , A4 = x 2 , X 5 = p 2 , 
A 6 = xp x +p x x, A 7 = y, A 8 = p y , A g = y 2 , Ai 0 = p 2 y , 
An = yp v +p y y, A12 = L z = xp y - yp x , A13 = xp y + yp x , 
A14 = xy and A15 = p x p y . The algebra exhibited by this 
set of operators is shown in Appendix A3. 

As in the previous examples, we first deal with the 
linear terms through the two-dimensional generalization 
of the transformation shown in Eq. (48) 

Ui = UitUixUiy, (169) 

where the t, x and y parts are given by 


V. TWO DIMENSIONAL CHARGED PARTICLE 
IN TIME-DEPENDENT ELECTRIC AND 
MAGNETIC FIELDS 

In this section we show that the general method pre¬ 
sented in Sec. II can be extended to obtain the evolu¬ 
tion operator corresponding to the Hamiltonian of a two- 
dimensional charged particle (—e) confined to a quadratic 
potential subject to an in-plane electric field and perpen¬ 
dicular magnetic field. The Hamiltonian of such a system 
is given by 

- 1 9 1 9 

H = — (p x + eA x ) + — {p y + eA y ) - ecf) 

2m 2m 

+ ^K(x 2 + y 2 )- (164) 


U u = exp 

I s « 




U lx = exp 

t n x it) x 

h 

exp 

i / \ * 
{t)Px 
h 

Ui y = exp 

y it) y 

exp 

— Xy (t) Py 


The corresponding transformation rules are 

UlPtUl = Pt + S AajIIaj XyTly 

+ X x Px H - H - ^yPyi 
UixUl = x + X x , 

Uiyill = y + X v , 

UiPxUl = Px Id-sc ? 

U\PyU\ = Py Ily- 


(170) 

(171) 

(172) 


(173) 

(174) 

(175) 

(176) 

(177) 
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Under this transformation the Floquet operator takes 
the form 




+ 


e 2 B 2 \ , 2 ^ 2 \ eB „, 

K + ^r) (x +y) + ^fv-vt‘*) 


— x 


' d dL 

dL' 

A 

' d dL 

dL' 

dtd X x 

dX x 

y 

dt dXy 

-1 

1 5>s 


+ Px 


' d dL 

dL ' 

+ Py 

' d dL 

dL ' 

. dt dti x 

dU x 

dtdtly 

-1 

5* 

C 


-Pt + L-S, (178) 


rotation of the form 
i 


U 2 = exp 


e ( t ) L z 


= exp 


-6 (t) ( Xf>y 


with transformation rules given by 


VPx ) , 
(185) 


&2PtUl =pt + 9 ( xp y - ypx) , (186) 

U 2 xU 2 = cos Ox — sin Oy, (187) 

fj 2 yU\ = sin#® + cos Oy, (188) 

UipM = cos 0p x — sin 0p v , (189) 

U 2 p y Ul = sin 0p x + cos 0p y . (190) 


Under U 2 the Floquet operator is transformed into 


yielding linear terms proportional to the Euler equations 
arising from the classical Lagrangian 




+ A^ILe + Aylly + eE x X x + eEyXy. (179) 


In order to eliminate the linear terms in x, y, p x and p y 
we demand that the Euler equations vanish 


d dL 

dL 

dt dX x 

dX x 

d dL 

dL 

dt dXy 

dXy 

d dL 

dL 

dtdtl x 
d dL 

a £ 
cS ^ 

dtdtly 

dBy 


- K 


e 2 B 2 
4 m 


—II y + U x — eE x — 0, (180) 
2 m 


- K + 


e 2 B 2 

4m 


Xy 


- - -IF,; + lly — eEy — 0, (181) 

-Ax - — - ^-Xy = 0, (182) 

m Zm 

-X y - ^ + ^-X x = 0, (183) 

m 2m 


and S = L. The Floquet operator becomes 



The third term is proportional to the z projection of the 
angular momentum L z = xp y — yp x which is the genera¬ 
tor of rotations around the z axis. We also notice that, 
besides the angular momentum L z , the first two terms 
given by the kinetic and potential energy are also invari¬ 
ant under rotations. Hence the next transformation is a 


Here it is important to stress that all the elements in the 
previous Floquet operator are invariant under rotations 
and therefore, the only extra element introduced by the 
transformation arises from the energy operator transfor¬ 
mation rule. By setting the restriction 

0 = eB/2m , (192) 

on the angle of rotation, the Floquet operator is trans¬ 
formed into the one of two uncoupled harmonic oscilla¬ 
tors 


0,0,(6-t.)0M = ±£ + l(K + ?£) S > 

In 1 / e 2 B 2 \ n 

+ 2^ p y + 2 [ K + ^r) y ~ Pu (193) 

with time-dependent parameters a = 1/m, 6 = 0 and 
c = K + e 2 B 2 /4m. 

It is clear how to proceed further: By repeating the 
procedure for the general quadratic Hamiltonian for the x 
and y harmonic oscillators. Even though in principle the 
two sets of transformations presented in Secs. IV A and 
IV B are equivalent, one is more effective than the other 
depending on the symmetries of the system. In cases, 
such as the charged particle subject to time-dependent 
magnetic field with varying mass, the set of transforma¬ 
tions of Sec. IV A yields closed and simple expressions 
for the transformation parameters whereas the transfor¬ 
mations of Sec. IV B give very complex ones. However, 
some other systems, as the charged particle with variable 
mass subject to constant magnetic field, are reduced us¬ 
ing a smaller number of simple parameters by means of 
the transformations of Sec. IV B. In the latter case, solv¬ 
ing Riccati differential equation is conveniently avoided 
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through the Arnold transformation. These two cases are 
presented as examples at the end of this section. 

Let us now reduce the Floquet operator (191) through 
the set of transformations from the first path (Sec. IV A). 
In this case, the evolution operator is given by 


& = ululutuluiul 


(194) 


where U\ and U 2 are given by Eqs. (169) and (185) re¬ 
spectively. The remaining transformations are general¬ 
izations of Eqs. (61), (67), (76) and (83) 


I / 3 = exp 


TffA (xpx + PxX) 


x exp 


^7 iVPy+PyV) 


t /4 = exp 
C /5 = exp 


—aAx 2 

2 h 


exp 


—a Ay 2 
2h y 


2h 


(j) ( xp x + p x x) 


x exp 


2h 


<t> iVPy + p v y) 


Uq = exp 


2W A 


exp 


—p?* 

27 A 


(195) 

(196) 

(197) 

(198) 


Since the x and y parts of the Floquet operator are sym¬ 
metric, the x and y part of these transformations have 
the same parameters and they may be obtained from the 
ordinary differential equations ( 66 ), (72), (81) and ( 88 ). 

If instead we follow the procedure from Sec. IV B the 
evolution operator is expressed as the product 


m = ululululululu\ 


(199) 


where, even though U\ and C /2 are newly given by (169)- 
(172) and (185), the remaining operators correspond to 
the generalizations of the transformations from the sec¬ 
ond path (61), (104), (113), (120) and (127) 


C / 3 = exp 


2h 


7 (xp x + p x x) 


x exp 


C /4 = exp 


T^liVPy +PyV) 


1 


2 h 


—<t> Ax + — p: 


x exp 


1 

2h 


1 


—( At/ + —p y 


C /5 = exp 
Uq = exp 


—aAx 2 

2 n 


exp 


—aAy 2 
2h y 


2h 


p ( xp x + p x x) 


x exp 


2h 


Pii)Py +PyV) 


C /7 = exp 


2n A 


exp 


—p V -A 

2U A 


( 200 ) 

( 201 ) 

( 202 ) 

(203) 

(204) 


The corresponding parameters may be obtained from the 
ordinary differential equations (103), (109), (118), (125) 
and (132). 

Having derived the explicit form of the evolution oper¬ 
ator (199) we obtain the Heisenberg picture position and 
momentum operators as 


' XH it) 


" X 


Xx 

m if) 

PxH it) 

= M 

y 

Px 

+ 

Xy 

-n x 

- PyH it) . 


-Py - 


- ~^y - 


(205) 


In this case, M is a 4 x 4 matrix that has the following 
form 


with 


GqqR G 9 pR 
GpqR Gpp R 


R = 


cos 8 — sin 8 
sin 8 cos 8 


(206) 


(207) 


the rotation matrix coming from the second transforma¬ 
tion. Notice that, as in previous examples, M and R 
have a symplectic form. The parameters A x , X y , n x and 
U y are calculated from the classical differential equations 
of motion (180)-(183), the rotation angle 8 is given by 
(192) and the coefficients G qqi G qp: G pq and G pp can 
be obtained by substituting the transformation parame¬ 
ters in (98)-(101) for the transformations in Sec. IV A or 
(136)-(139) for the transformations in Sec. IVB. 

For the first path, the Green function is calculated by 
placing the transformations (169), (185) and (195)-(198) 
in Eq. (15) and using the explicit form of the propagators 
presented in Appendix B 


G (x,y,t;x' ,y' ,0) = J dxidyidx 2 dy 2 ... dx 5 dy 5 


>x,y 

x (x 2l y 2 

x exp 


ul 


xi,yi) ( 27 , 2/1 


Ul 


£ 2 , 2/2 


U. 


£3,2/3 


£ 5 , 2/5 


Ul 


x',y' 


2nh/3 


e ^ 7 exp I — * 


(£ - A x ) - ( 2 / - A y ) 


x exp < i 


.Ae" 27 / e~ 2 ^ 


2h 


p 


■ — a 


(x— A x f + (y- A V Y 


x exp < 1 


A 

2hj3 L 


,/\2 


(£') +( 2/0 


Ae _,i_7 

x exp •( —i ——— \(x — A x ) (x' cos 8 — y' sind) 

l 

+ ( 2 / — A y ) i x ' sin 8 + y' cosd)]| . (208) 

Similarly, the second path’s propagator is calculated by 
;athering the explicit form of the seven transformation 
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propagators given in the Appendix B and performing the 
integral in Eq. (15). This procedure yields 

G (x, y, t; x ', y', 0) = J dxidyidx 2 dy 2 ■ ■ ■ dx 6 dy 6 

u\ 


x,y 


x (*2,2/2 


* 1 , 2 / 1 ) (* 1 , 2/1 


m 


* 2 , 2/2 


u\ 


*3,2/3 


*6,2/6 


U] 


*', 2 /' 


A exp (—ip — 7 — iS/h) 


x exp |zu ^(aj - A x f + (y - A y ) 2 J | 


x exp 


iq ( x' cos 9 — y' sin 6 — 2/- ) (x — \ x ) 
qh) 


IT, 


+iq ( x' sin 9 + y' cos 9 - Til ) (y — \ y ) 


, (209) 


where the functions w, u, q and l are given by Eqs. (141)- 
(144). 


A. Charged particle in a time-varying magnetic 
field 


Here we consider a charged particle subject to a mag¬ 
netic field of the form B = Bo smut. Since there are 
no linear terms A^, = A y = Hr = I \ y = S = 0. By in¬ 
troducing this form of the magnetic field in Hamiltonian 
(164) and applying the first two transformations (169) 
and (185) we are left with the Hamiltonian of two un¬ 
coupled harmonic oscillators 

&£>, (h - #,) om = 2 '-ti+ YJt * 2 


+ 


2771 


„ 2 1 e 2 Sg 

v A -- 

Py 2 4m 


sin 2 ut y 2 — p t , (210) 


where we identify a = 1/m, b = 0 and c = 

e 2 Bo sin 2 ut/Am = e 2 Bo (1 — cos 2ut) /8m. Given that 
in this example the mass does not depend on time and 
consequently a =cnt. then 7 = 0. Placing a, b and c 
in the ordinary differential equations (72), (94), (95) and 
(192) we obtain the following solutions 


( Suj 2 ’ 16tZ 2 ’ w< ) uC' (ut) 

f ui 2 Ul 2 
8t^ 2 ’ 16u 2 1 L 


c 


“ =- - - <211) 

</ = ln c [J^, Tfa2 ,UJt ) = ln [CM)], ( 212 ) 

13=1 , ,„ a - ~=T(ut), (213) 

= — sm-—. (214) 

UJ Z 

where u c = eBo /m and C and C' are the cosine Mathieu 
function and its derivative which comply with C (0) = 1 
and C' (0) = 0. 


2 vt 


U c 


The Heisenberg picture position and momentum oper¬ 
ators are obtained by replacing the previous parameters 
in (98)-(101) and (205) 


xr = C (ut) 

T (ut) 


U c .9 Ut 
x cos | — sm — 
u 2 


c . 2 ait 

V sm ( — sm — 
UJ z 


C(ut) 


U c 2 ait 
p x cos I — sm — 
u 2 


c . 2 ait 
-Py sm ( — sm — 


( 215 ) 


y H = C (ut) 

T (ut) 


ut 


x sm — sm — + y cos — sm — 

u 2 J \ u 2 ) 


ut\ 


m 


C (ut) 


U c . 2 Ult 
p x sm [ — sm — 
u 2 


. . . U c . 2 Ut 

+Pv cos — sm — 
' u 2 


PxH = 


1 


+ uC' (ut) T (ut) 


C (ut) 

„ . f U c .2 Ut 

-p v sm — sm — 
y UJ Z 

^ • f ^c 2 ^^ 

~y sm —sm — 

V to 2 


U c , 2 u t 
p x cos ( — sm — 

UJ Z 


A- muC' (ut) 


U c . ' 

x cos | — sm 2 
u 


PyH = 


C (ut) 


+ uC' (ut) T (ut) 


c . 2 

p x sm [ — sm — 
UJ z 


. ~ , c . 2 ^ 

+p y cos I — sm — 


UJ 


~ c . 2 

+y cos — sm — 
cj 2 


• mujC' (ujt) 


x sm I — sm 

UJ 


(216) 


ut 

T 

(217) 


ut 

T 

(218) 


The propagator is calculated by introducing the ex¬ 
plicit forms of the transformation parameters in (208) 
giving 


G(x,y,t-,x',y’,0) = 
x exp 


2nUT (ut) C (ut) 

1 uC' (ut) 


2tyU \T (ut) C 2 (ut) C (ut) 


U'? + (y'Y 


(* 2 + y 2 ) 


x exp < 1 


x exp < —1 


2TiT (ut) 

m I". , (u c . 2 

(xx + yy ) cos — sm — 

L \u 2 J 


UT (ut) C (ut) 


. 1 j. 1 (Uc 2 ait 

+ (yx — xy ) sm ( — sm — 


(219) 


B. Charged particle in time-dependent electric 
fields 


The example treated in this section may be of use in 
modelling single electron quantum dots [44], or magneto 
transport in semiconductors under the influence of an in¬ 
cident radiation [4], In the latter application, the degree 
of circular polarization plays an important role that may 
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be elucidated through the model presented in this sec¬ 
tion. In order to allow the possibility of studying the 
effects of polarized incident light we introduce the fol¬ 
lowing form of the electric field 


E x = E 0x + E lx sin (ut ), (220) 

Ey = Egy + Ely Sill (tilt + ^) , (221) 


in Hamiltonian (168) where Eo x and E^y may be consid¬ 
ered bias electric fields and ( controls the degree of cir¬ 
cular polarization of the incident radiation with electric 
field components E\ x and E\ y . The first transformation 
to perform is (169) that yields the classical differential 
equations of motion (180)-(183). The solution to these 
equations is obtained after a lengthy calculation 


Aa; — 


4eE 0x lQeEi v ujui c cos ( cos (tu) 


sin(fu;) 


T±m 

16eEi x ui 2 


r 4 


m 


+ sin 


Sill 


tu c 

~2 


T 4 m 

4Fi eEi x 16eEi y uu} c sin £ 
T 4 m T 4 m 

4eE 0x u c 32eEi y uj 3 cos ( 


+ cos 


T 2 _mQ r 4 mQ 

8eEi y uu 2 cos ( 8eEi y uflcos(' 

T 4 mf2 
(tu c 


m 

— cos 


V 2 


T 4 m 

/ 4eEo x 16eEi y uu c cos ( 


- cos 


Sill 


tu c 

~2 


\ Tim 

4eEo y 


T 4 m 

16eEi x uu c 


+ 


rim T 4 m 

16eEiy0j 2 sin ( 4rieE\ y sin 
T 4 m 


Sill 


T 4 m 

/ 4eEoyto c 32eEi x uj 3 

r 4 mfl 


V ri mfl 


m\ (tu c 
~2 ) C ° S V~2" 

8r+eEi x uj l6eEi y u 2 u c sin £ 

^ r 4 m n r 4 mfi 

4rieEi y u>c sin£ 
T 4 mfi 


( 222 ) 


^v ~ p2 


4eEo y ^ ^ f 16eE\ x uu c 


+ cos (tu) I — 


m V T 4 m 

l6eEi y u; 2 sin £ 4rieE\ y sin £ 


T 4 m 


T 4 m 
^2 


/ 16eEiyUj 2 cos i 4r 2 _eEi y cosC\ 

^ V r 4 ™ f*m ) 

4eEo x u c 32eE\ y u 3 cos ( 


'tQ\ (tu, 
+ S inl T lco s l T 


rimfl 


r 4 mfi 


8eEi y ojuj 2 cos £ 8eEi y uQ cos£ 

+ r 4 mfl T 4 m 

tn\ . (tu c \ f4 cEq x 16eEi y uu c cos( 

2 ) S \ 2 j \ rim r 4 m 

(tQ\ (tu c \ (4eEo v 16eEi x uu c 

\2 J V 2 / V r - m T m 

IQeEiyU 2 sin ( 4r 2 _eEi y sin ( 


T 4 m 


+ 


T 4 m 


. ftn\ . (tu c \ (4eE 0y u c 32eEi x u> 3 

\ 2 J \ 2 J \ rimfl T 4 mf2 

8r 2 ^_eEi x uj 16eEi y u 2 u c sinC, 


r 4 mVl T 4 mQ 

4TieEi y ijj c sin£ 

r 4 mfl 


, (223) 


n* = 


2cE(ly(jJ C 


ri 


+ sin(fw) 


( 8eEi y u 2 u c cos £ | 2TieEi y u c cos £ 

I p4 p4 

16e£i x u; 3 4r 2 ,eEi x u 


+ cos (tui) 

8eEi y u 2 u c sin £ 2rieE\ y u c sin£ 


+ 


r 4 


r 4 


f tO,\ . f tu c \ ( 2eE 0x u c IQeEiybj 3 cosC 

+ cos (Y) sm ( _ 2“) ( rl-fi- 

4eEi y ujQ 2 cos £ 4eE\ y uu 2 cos £ 


T 4 


r 4 


'tfl\ (tu c \ (2eEo x n 8eEi v uQu c cos( 

sm| Y lc° S l—j [-=-p- 


'tn\ . ftu, 

sm ( — j sin ( — 


2eEo v n 8eEi x uQu> c 


ri r 4 

8eSi y o; 2 HsinC 2T?_eTiiyf2sinC s 


T 4 


r 4 


'tfl\ (tu c \ ( 2eEo v u c 16eE\ x uj 3 
+ co B | —Icosf — 1 t -+ 


4r 2 Sr eEi x u 8eE\yU 2 u c sin £ 

p4 p4 

ZTieEiyUc sinC 

f 4 


. (224) 
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i, ^eEflxUJc 

ll y ~ p 2 


+ sin(£w) ( — 


8eEi x u 2 iu c 2T 2 _eEi x uj c 


+ 


p4 p4 

16eEi y oj 3 sin£ 4eEi y ujfl 2 sin ( 4eEi y ujw 2 sin£ 

p4 p4 p4 

4T+ eEiyto cos ( 16eEi y uj 3 cos £ 


+ cos (tuj) 


p4 


p4 


{ttt\ (tu) c \ f 2eEo x uj c 16eEi v Lo 3 cos C 
+ co^ T jco ^—j [~^—+ -JI- 

AeEiyUifl 2 cos ( AeEiyUJW 2 cos £' 


p4 


p4 


. . (tuj c \ (2eEo x fl 8eEi v ujfluj c cos( 

+ smi Y>‘ m {—){—_ -F3- 

' tfl\ (tco c \ (2eEo v Q 8eE\ x uiflu> c 

+ sin I — cos 1 11 


2 J \ 2 J \ T 2 _ r 4 

8eEi y uj 2 fl sin £ 2T 2 _eEi y flsin( 


p4 


p4 


(tfl\ . ( tu c \ ( 2eE 0v u> c lQeEi x uj 3 

+ co s | T U„l— j I-— + ^ 3 — 


AT^eEixU 8eEi y u} 2 Lu c sin £ 

p4 p4 

2 T/. eE\y uj c sin£ 
T 4 


. (225) 


where ui c = eB/m is the cyclotron frequency, 
il 2 = (4 K + mix 2 ) /to, T 4 = 16w 4 + (fl 2 — oj 2 ) 2 — 
8 w 2 (fl + w c ) 2 , T 2 + = fl 2 + co 2 and T?_ = Q 2 — ix 2 . 

Integrating Eq. (192) yields the rotation angle given 
by 


„ eB ix c 
6 = -z—t = -7rt. 
2m 2 


(226) 


x H 


- UU£> 


2 (fi \ 
H-— sin ( — t 

toO \ 2 J 
4~X x , 


/ ix c \ 
Px cos tj 


(231) 


VH 

+ Xy, 
Px = CO 


(fl \ 

[A • C w c , 

A , - C w c 1 

\1 


ism t 

J + 2 /cos (^y t 

)J 


2 . 

(fl A 

r. . /w c , 

A , - (u c 

\1 

+ M sm 

V 2 f ) 

PiSin (^y t 

j +P y COS ^y I 

)\ 


(232) 


(fl A 

r. (u c 


\1 

U*) 

Px COS (^y t 

j -Pj/Sin^y t 

)\ 


mfl ( fl \ 

sm U *) 

-fix, 


£cos(y tj -ysin(y 


(233) 


p v = c 

mfl 

2 

~bly 


(fl A 

r. . /We . 


\1 

U 4 ) 

Px sm ^y £ 

) +Py cos^y f 

)J 


sin ^ tj x sin (y tj+y cos (y tj 

-ry (23 

In the calculation of the propagator, we only consid 
the transformations t/i, U 2 and C /4 since they have no 
vanishing parameters X x , X y , fl x and 11^ and A. I 
inserting these transformations in (15), the propagat 
takes the form 


G {x, y, t, x , y , 0 ) , . / q , 

4nh sin (j t) 

' ” 1 r \ , 

T~(y~^y 


-iS/h 


x exp 


■ fix / , S 

- 1 — (X - Xx) 


exp 


,n„ 


Xy 


As mentioned above, this equations transform the Flo- 
quet operator into the form of two uncoupled harmonic 
oscillators with parameters a = 1 /m, 6 = 0 , c = 
K + m(w c /2) 2 , d = 0, e = 0 and g = 0. 

In this example we follow the sequence of transforma¬ 
tions from the second path presented in Sec. IV B which 
rapidly reduce the Floquet operator avoiding the solution 
Riccati equation in contrast to the procedure presented 
in Sec. IV A that yields a larger number of parameters. 
The remaining parameters 7 , </>, a, ip and fd are obtained 
from Eqs. (103), (109), (118), (125) and (132). These 
equations yield 


7 = 0, 

(227) 

mfl 

A = —’ 

(228) 

, O 
^=2 

(229) 

a = ip = /3 = 0. 

(230) 


x exp < 1 - 


mfl 


{%'? + W? 


, ft 

cos — t 


46, sin (y t) L 
+ (x — X x ) + (y — X y ) 

- 2x' (x - X x ) cos (y tj - 2x' (y - X y ) sin (y tj 
+ 2y' (x-X x ) sin (y tj - 2 y' (y - X v ) cos (y tj |. 


(235) 


VI. CONCLUSIONS 

The Lie algebraic technics relay on the existence of a 
set of generators that forms a closed algebra. If a given 
Hamiltonian can be expressed as a linear combination 
of these generators, the overall structure of its evolution 
operator is known and takes the form of Eq. (5). 
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We have applied the Lie algebraic approach to obtain 
the evolution operator of the general harmonic oscillator 
and the charged particle in time-dependent electric and 
magnetic fields. The sets of operators that form closed 
Lie algebras characterized by their structure constants 
where established in each case. Some particular exam¬ 
ples of these two Hamiltonians were examined in detail. 
The free particle in the presence of an external driving 
force was used to introduce the Lie algebraic approach. 
Analytical expressions for the evolution operator were 
provided for the potential of a radio frequency ion trap 
as well as for a forced harmonic oscillator with varying 
mass. The charged particle’s evolution was studied under 
two sets of different conditions. First we treated the case 
of a sinusoidally varying magnetic field and second, we 
calculated explicit expressions for the evolution operator 
and propagator of a particle in constant magnetic field 
and time-dependent electric field. 

The here obtained results may be used to tackle di¬ 
verse problems as squeezed states, radio frequency traps, 
and electronic transport in two-dimensional lateral het¬ 
erostructures under diverse conditions of light excita¬ 
tion (linear polarization, circular polarization, etc.). The 
methods developed so far can be extended to other types 
of Hamiltonians, for example, a charged particle subject 
to time-dependent electric and magnetic fields in a asym¬ 
metric confining parabolic potential. 

We have observed that the Lie algebraic approach is a 
powerful method that can be used to obtain the evolution 
operator and propagator of a great variety of Hamiltoni¬ 
ans. It reduces the the difficulty of solving Schrodinger 
partial differential equation into obtaining the solution of 
a system of coupled ordinary equations for the transfor¬ 
mation parameters. Two possible shortcomings of this 
method are that in general the solution of quite com¬ 
plex ordinary differential equations is needed, and more 
important, it requires a finite dimension set of operators 
that form a closed algebra which in many cases is difficult 
to identify. 
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Appendix A: Generators of the Lie algebra 

In this appendix we list the generators, and their cor¬ 
responding Lie algebras characterized by the structure 
constants for the three examples treated in this paper: 
the linear potential, the general quadratic Hamiltonian 
and the Hamiltonian of a charged particle subject to elec¬ 
tromagnetic fields. 


1. Generators for the linear potential Hamiltonian 

The set of Hermitian operators that form the closed 
Lie algebra for the Hamiltonian (18) is given by 


= 1, 

(Al) 

= X, 

(A2) 

= p, 

(A3) 

=f , 

(A4) 


where 1 is the identity operator. The commutation re¬ 
lations arising from all the possible combinations of the 
previous generators yield 


Ai 

A2 

= M 

= 0, 

(AS) 

Ai 

A3 

= M 

= 0, 

(A6) 

Ai 

a 4 

= M 

= 0, 

(A7) 

A2 

A3 

= [x,p] = 

ihl = ihXi, 

(AS) 

A2 

a 4 

= [x,P 2 ] 

= iTi2p = ih 2A3, 

(A9) 

A3 

a 4 

= [p,p 2 ] 

= 0, 

(A10) 


therefore the structure constants are 02,3,1 = 1, 02,4,3 = 2 
all others begin zero [see Eq. (3)]. 


2. Generators for the general quadratic 
Hamiltonian 

The structure of the quadratic Hamiltonian (47) sug¬ 
gests that the closed algebra is given by the set of oper¬ 
ators 


Ai — i, 

(All) 

A 2 = x, 

(A12) 

A3 =P, 

(A13) 

A 4 = i 2 , 

(A14) 

As =P 2 , 

(A15) 

Xq = xp + px. 

(A16) 
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Indeed, the commutation relations for these operators 
yield a closed algebra given by 


Ai, A2 

= 

Ai, A3 

= 


= 

Ai, A6 

= 


— Ai,A 4 — Ai,A, 


A2, A3 
A2, A4 
A 2 , A 5 
A 2 , Ag 
A3, A4 
A3, A5 
A3, Ag 
A4, A5 
A4, Ag 
A5, Ag 


i, AiJ = 0, (A17) 

[x, p] = ihl = ihXi, (A18) 

[i, x 2 ] = 0, (A19) 

\x,p 2 ] = ih2p = ih2X 3l (A20) 

[i, xp 4 px] = ih2x = ih2\2, (A21) 


[p,x 2 

[p,f] - 

\_P, xp 4 

[x\f] 

[i 2 , xp 
[p 2 ,xp 


= —ih2x = —ih2X.2, (A 22 ) 

= 0, (A23) 

- px] = —ih2p = —ih2X 3 , (A24) 

= ih2 (xp 4 px) = ih2Xe, (A25) 
4 px] = ihAx 2 = ihAX±, (A26) 
4 px] = —ihAp 2 = —iMAs.(A27) 


The structure constants are 02 , 3,1 = 1, 02 , 5,3 = 2, 02 , 6,2 = 
2 , 03 , 4,2 = — 2 , 03 , 6,3 = — 2 , 04 , 5,6 = 2 , 04 , 6,4 = 4 and 
05 , 6,5 = ~4, all others being zero. 


3. Generators for the Hamiltonian of a charged 
particle in electromagnetic fields 

The structure shown by the Hamiltonian of a charged 
particle in electromagnetic fields (168) suggests that the 
set of generators that yields the corresponding closed Lie 
algebra should be composed of the identity operator, the 
generators listed in the previous section (A 2 to Xq) for the 
x and y parts of (168) plus the angular momentum L z = 
xpy — ypx- However, the algebra formed by this set needs 
three more operators in order to be closed. Thereby, 
the whole set of operators is given by Ai = 1 , A 2 = x, 
A3 = Px, A 4 = x 2 , A 5 = pi, A 6 = xp x 4 p x x, X 7 = y, 
As =P V , Ag = y 2 , A10 =p 2 , A n = yp v +p v y , A 12 = L z = 

xp y ~ VPx, A 13 = xp y 4 yp x , A 14 = xy and A 15 = p x p y . 
Here is a summary of all the connnutors arising from the 
generators listed above, for convenient reference: First, 
all the generators commute with the identity operator 


Ai, A2 


Ai, A3 


Ai,A 4 — Ai,A, 


Ai, As 


Ai, A 7 
Ai, An 


Ai, Ag 
Ai, A12 


Al,Ag = 

Al,Ai3 = 
= Ai,Ai5 


Ai, A10 
Ai, A14 
= 0. (A28) 


Second, any generator from the x part commutes with 
any generators from the y part 


A2, A 7 



A 

2, As 

= 

A: 

, Ag 

- 

A2 

, A10 


= 

^2 

, An 


A3, 

A7 

= 

A 3 , 

As 

= 

A3, A9 

= 

^3 

A 

10 

= 

A3, An 


A4 

A7 

= 

A4, As 

= 

a 4 

A 

9 

= 

A4, A 

10 

= 

A4, An 

= 

A5, a 7 

= 

A5 

A 

8 

= 

A5, A9 


^5, A 

10 

= 

A5, An 

= 

^6 

, a 7 

= 

A6, As 

= 

A6, Ag 


^6, A10 


A6, An 


= 0. (A29) 


Third, the generators belonging to the x coordinate must 
follow the commutation rules (A18)-(A27), therefore 


A2, A3 
A2, A4 
A2, A5 
A2, A6 
A3, A4 
A3, A5 
A3, A6 

A4, A5 


= [£,Pa;] = iTi 1 = ihXi, 

= [x,x 2 ] = 0 , 

= [x,p 2 ] = ih2p x = ih2X 3 , 


= [x, xp x 4 p x x] = ih2x = 1 T 12 X 2 , 


= [pa,,! 2 ] = —ih2x = —ih 2 X 2 , 
= [Px,pI] = 0 , 

= \Px,xp x 4 p x x] = -ih2p x 
= —ih 2 A 3 , 

= [x 2 ,p 2 ] = ih2 (xp x 4 PxX) 


(A30) 

(A31) 

(A32) 

(A33) 

(A34) 

(A35) 

(A36) 


A4, A 6 
A5, A6 


= ih2X 6 , (A37) 

= [x 2 ,xp x 4 p x x] = ihAx 2 = ihA A 4 , (A38) 

= [p 2 x, xpx + Pxx] = —ihAp 2 

= -ihAX 5 . (A39) 


Similarly for the y coordinate we have 


a 7 

00 

1 - 1 

v> 

, Ag 

a 7 , 

A10 

a 7 , 

An 

As 

, A 9 

00 

A10 


[y,p v ] = ih 1 = ihX-L, 

(A40) 

to 

II 

0 

(A41) 

[y,P 2 y\ = ih2p y = ih2X 3l 

(A42) 

[y,yp y +p v y] = ih2y = ih2X 7 , 

(A43) 

[py,f] = -ih2y = -ih2Xr, 

(A44) 

[p V ,P 2 y] = 0, 

(A45) 
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As, An 

= [Py,yPy+P y y} = -iTttpy 


A 6 , A 12 


= —ih2X 8 , 

(A46) 


i- 1 

to 

V> 

o 

= [y 2 ,p 2 v \ = ihz (yp v + p v y) 


A 6 , A 13 


= ih2Xn , 

(A47) 


1 1 
V> 

CO 

= [ y 2 ,yp y +Pyii] = iMy 2 = 

iMA9,(A48) 

A 6 , A 14 

Aio, An 

= \p 2 y ,ypy+p v y\ = -iMpl 




= —iMAio- 

(A49) 

A6j Ai5 


The remaining operators need to be calculated indepen¬ 
dently 


A 2 , A 12 

= 

x, L z 

= —ihy = —ihX 7 , 

(A50) 

A 7 , A 12 

A 2 , Al3 

= 

x, xp y + yp x \ = ihy = ihX 7 , 

(A51) 

A 7 , A 13 

A 2 , A 14 

= 

x,xy] 

= 0 , 

(A52) 

A 7 , A 14 

A 2 , A 15 

= 

X,P X Py] = iliPy = ihXs, 

(A53) 

A 7 , A 15 


A 3 , A 12 

= 

Px,L z 

= —ihpy = —ihXs, (A54) 

As, A 12 

•A 3 , A 13 

= \p xi xp y + yp x ] = -ihp y = -ihx 8l ( A55) 

As, A 13 

A 3 , Ai 4 

II 

= —ihy = —ihXy, (A56) 

As, A 14 

A 3 , A 15 

= [PxiPxPy} = 0, (A57) 

As, A 15 


XPx + PxX, L z 


—ih2 ( xp v + yp x ) = -ih2X 33 , 

(A 66 ) 

[xp x + p x x, xpy + yp x \ 

—ih2L z = —ih2X\2, 

(A67) 

[xp x +p x x,xy} 

—ih2xy = —ih 2 Xi 4 , 

(A 68 ) 

[Xpx +PxPx,PxPy] 

%h2p x Py - lh2X\$, 

(A69) 

y, L z = ihx = ihX 2 , 

(A70) 

[y, xpy + ypx] = ihx = ihx 2 , 

(A71) 

5cg> 

II 

O 

(A72) 

[y,PxP y ] — ihp x = ihx 3 , 

(A73) 


p y ,L z = ihp x = ihX 3 , 

[py, xp y + yp x \ = -ihp x = 

[p y ,xy\ = —ihx = —ih\ 2 , 
[Pyi PxPy] — 0, 


(A74) 
-ih\ 3 ,(A75) 
(A76) 
(A77) 


A4, A12 

= x 2 ,L z = — ih2xy — — ih2X\4, (A 58 ) 

Ag, A12 

A4, A13 

= [i 2 , xp y + yp x \ = ih2xy = i?i2Ai4j(A59) 

Ag, A13 

A4, A14 

= [x 2 ,xy] = 0 , (A 60 ) 

Ag, Ai4 

A4, A15 

= [x 2 ,p x p y \ = ih2xp v 

Ag, A15 


= ih ^Ai 2 + Ai 3 ^ , (A61) 



ih2xy = ih2Xi 3 , 


(A78) 


= [y 2 , xp y + yp x ] = iU2xy = i7j,2Ai 4 ,(A79) 
= [f,xy\= 0, (A80) 

= [y 2 ,PxP v \ = ihypx 


= iTi[\ 


13 



(A81) 


A5, A12 

= 

' ' 

H to 

= —ih 2 p x p y = — ih 2 X 33l (A62) 

A10, A12 

•A5, A13 

= 

p 2 x , xpy + yp x \ 

v>' 

0 

v> 

CO 


= —ih2p x p y = —ih2\i 3 , (A63) 


An, A 


14 


= [pl,xy\ = -ihZypx 




= ih 




(A64) 


An, A 


15 


= [pl,PxP y \ = 0, 


(A65) 





ih2p x p y = ih2\i 3 , (A82) 


\p 2 y,xp y + yp x ] 

—iTi2p x p y = —ih2\i 3 , (A83) 

[p 2 y ,xy\ 

—ih2xp y = —ih ^Ai 2 + Ai 3 ^ , (A84) 
[p 2 v ,PxPy\= 0, (A85) 
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An, A 12 

An, A13 

An, A14 

An, A15 


= yPy+PyViL^ 

= ih2 (xpy + yp x ) = ifi.2Ai 3 , (A86) 
= [yPy +Pyipxp v + PPx\ 

= ih2L z = , (A87) 

= [ yp v + p v y , xy\ 

= —ifi2xy = —i?i2Ai4, (A88) 

= [ypy +p y y,PxP v \ 

= ih2p x p y = ih2\is, (A89) 


Al2, A13 


Al2, A14 


Al2, A15 


= L z , xp y + yp x 
= ih (VPy + PyV - XPx - PxX) 

= ih ^An — XqJ , (A90) 

= L z ,xy 

= ih (y 2 — x 2 ) = ih ^Ag — A 4 ) , (A91) 

= L z , P X Py 

= ih (•py - pl) = ih ^Aio - A 5 ) , (A92) 



A13 

, A14 

II 

+ yPx,xy\ 






= —ih 

( x 2 + y 2 ) = —ih 

^4 + Ag^j ,(A93) 


A13 

, A15 

= [xpy 

+ VPx,PxPy} 






= ih (/ 

>1 + P 2 y) = ih (X t 

+ Am) , (A94) 


A14 

, A15 

= [xy, 

PxPy\ 






= ih— 
2 

(■XPx + PxX + yp v + Py 7) 





= ih— 
2 

^Ag + An^j . 

(A95) 

The 

structure constants inferred from 

the 

15!2! / 

(15 

— 2)!2! = 105 commutors 

in Eqs. (A28)- 

(A28) 

are 

£2,3 

1 — 1) c 2,5,3 — 2, £2,6,2 

= 2, £3,4,2 = 

-2, 

c 3,6,3 

= — 

2, c 4 

5,6 — 2, 

£4,6,4 = 4, £5,6,5 

= —4, £7,8,1 

= 1, 

<-7,10,8 

= 

2, £7,11,7 = 

2) C8,9,7 = — 

2, £8,11,8 = 

-2, 

£9,10,11 — 

2, £9,11,9 = 

4) £10,11,10 = - 

-4, £2,12,7 = 

-1, 

C2,13,7 

= 

1) C2,15,8 = 

1) £3,12,8 = - 

1) C3,13,8 = 

-1, 

£3,14,7 

- 

-1, 

C4,12,14 

= —2, £4,13,14 = 

= 2, C4,15,12 

= 1, 

C4,15,13 = 

1, C5,12,15 - 

—2, £5,13,15 = 

—2, £5,14,12 

= 1, 

C5,14,13 = 

-1, 

£6,12,13 - 

= —2, C6,13,12 = 

-2, £6,14,14 = 

—2, 

C6,15,15 = 

2, 

£7,12,2 

= 1) £7,13,2 = 

1) £7,15,3 

= 1, 

C8,12,3 

= 

1) £8,13,3 — 

— 1) C8,14,2 = 

-1) C9,12,14 

= 2, 

C9,13,14 — 

2, Cg,i5,i2 = 

= —1) £9,15,13 = 

1) Cl0,12,15 

= 2, 

Cl0,13, 

15 

= 

-2, £10 

14,12 = —1) 

£10,14,13 = 

-1, 

£ll,12, 

13 = 

2, £11,13,12 

= 2, £11,14,14 = 

— 2, £11,15,15 

= 2, 

£12,13, 

6 = 

-1, 

Cl2,13,ll 

= 1, Cl2,14,4 = 

— 1) £12,14,9 

= 1, 

£12,15, 

5 = 

-1, 

Cl2,15,10 

= 1) Cl3,14,4 = 

_ 1) £13,14,9 = 

-1, 


Cl3,15,5 — 1, Cl3,15,10 — 1, Ci4 I 15 i 6 — 1/2, 014 , 15,11 — 1/2, 
all others being zero. 


Appendix B: Unitary tranformations 

The following sections are devoted to presenting the 
unitary transformations generated by the operators in 
the previous appendices used along the paper to reduce 
the different Floquet operators. The general form of 
each transformation is accompanied by the transforma¬ 
tion rules, i. e. the explicit forms of UptU\ UxU 1 and 
UpU'. The transformation’s Green functions are also 
presented in the sections to follow. 

Since it is widely used in the following appendices we 
enunciate the next commutation relation. If the commu- 
tor 


A, B =C, 


(Bl) 


commutes with the operators A and B , i. e. 


A. C 

then it follows that 

A,F(B) 


B,C 


A,B 


= 0 . 


dF{B) 
dB ’ 


(B2) 


(B3) 


provided that F is an analytical function. 


1. Unitary transformation generated by x and p. 
Shift in space and momentum 


This transformation shifts the space and momentum 
operators by time-dependent functions. It is generated 
by Ai, A 2 and A 3 in Eqs. (All), (A12) and (A13) as 
follows 


U = UtUxU P = exp 



x exp 


-II(f)a 


exp 


-A (t)p 


(B4) 


The transformation rules for the space, momentum and 
energy operators can be worked out by inserting commu- 
tors 


UxU" 1 = x + U P x,Ul 
UpU'=p + U x [p,U' x 
Up t U' = p t + Ut \pt, u' t 


(B5) 

(B 6 ) 

(B7) 
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and using relation (B3) as follows 


3. Transformation generated by x 2 


UxU^ = x + U P [x,p] 


dU$ 

dp 


— x + A , 


UpW =p + u x [p,x] -ul=p-n , 

Upt& =Pt + Ut\pt,t]^- t U} 

= p t + s — An + rix + a p . 


(B 8 ) 

(B9) 

(BIO) 


The propagator for this transformation is given by 
(x |[/t| x') = exp 


iS 

h 


exp 


m . 

-~h x 


x S (x — x' — A), (Bll) 
where S is the Dirac delta distribution. 

2. Transformation generated by xp + px. Dilation 

Dilations are generated by A 6 in Eq. (A16). The ex¬ 
plicit form of this transformation is given by 


U = exp 


2h 


7 (t) (xp + px) 


(B12) 


In order to get the transformation rules for the position 
and momentum operators we define 


xfr) = UxU\ 
P( 1 ) = UpU\ 


(B13) 

(B14) 


and compute the derivatives with respect to the trans¬ 
formation parameter 7 

^~X(l) = -^U [xp+px, x\U^ =X , (B15) 

= xrU [xp+px, p]U^ = -P . (B16) 

The solution to this pair of differential equations together 
with initial conditions 2 f( 0 ) = x and P( 0 ) = p yields the 
standard transformation rules for dilations 


X = x e 7 , 
P = p e -7 . 


(B17) 

(BIB) 


The transformation rule for the energy operator is easily 
calculated by inserting a commutor 


Up t U^ =p t + U p t ,l 7 t 

and using relation (B3) as follows 


(B19) 


d 


Up t tf =p t +ih U^-W =Pt + \ (xp + px) . (B20) 

The corresponding propagator is given by 

(x |[/ 1 ’| x') = e~i 5 (e _7 x — x') . (B21) 


Here we analyze the transformations generated by A4 
in Eq. (A14). 


U = exp 


Y(t) 


Ax 2 


2 h 


(B22) 


Since it depends explicitly on the position operator, the 
position operator itself remain unaltered under its action 


UxU 1 = x . 


(B23) 


The momentum and energy transformation rules are 
easily obtained by inserting a commutor and using rela¬ 
tion (B3) as follows 


UpU 1 = p + U\p , x] 


dW 

dx 


= p — a Ax 


(B24) 


/s A. F) * CV /\ 

Up t W =pt + ihU—W =p t + —x 2 . (B25) 

The propagator associated to this transformation is easily 
calculated by using (B24) 

(x |[/ i '| x') = exp S(x — x'). (B26) 


4. Transformation generated by p 2 

This transformation is generated by A 5 in Eq. (A15) 
given by 


U = exp 




(B27) 


Under the action of U, p remains unaltered since these 
two operators trivially commute 


UpW = p . 


(B28) 


The energy and position transformations rules are eas¬ 
ily worked out by inserting a commutor and using the 
relation (B3) as follows 


- ' - dW 

UxU^ = x + U\x,p]-—- = x + (3p , 

op 

UptU f = p t + in u^-w =p t + . 

at 2 

This transformation’s propagator is given by 

1 


(x |[/t| x') = 


^2nhp 


exp 


/\ 2 


2fi,/3 


(x — x') 


(B29) 

(B30) 

(B31) 
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5. Transformation generated by Ax 2 + p 2 /A. 
Arnold transformation 


The Arnold transformation is generated by a linear 
combination of A 4 and A 5 in Eqs. (A14) and (A15). This 
transformation’s explicit form is given by 


U ( t ) = exp 


JLWa * 2 *^ 2 


(B32) 


where A is a constant that yields a unit-less <j> transfor¬ 
mation parameter. 

The transformation rule for the energy operator is 
readily calculated by inserting a commutor as follows 


UptW =pt + u \p t , Eft] =p t + u [put] 


< 917 + 

dt 


= Pt + 


p 2 + Ax' 


(B33) 


In order to obtain the transformation rules for the posi¬ 
tion and momentum operators we define 


X{cj>) = Uxtfl, 
P(0) = UpU\ 


(B34) 

(B35) 


and compute the derivatives with respect to the trans¬ 
formation parameter 




x, Ax 2 + —p 2 U + 






p, Ax 2 

= ax(4>). 


A 1 


U + 


(B36) 


(B37) 


The solution to this system of differential equations to¬ 
gether with the boundary conditions X (0) = x and 
P (0) = p yields the transformation rules 

1 


UxW = x cos cj) + -^p sin (j>, 


(B38) 

Up& = pcoscf) — Axsin</>. (B39) 

The Arnold’s transformation propagator is given by 


(x | IT 1 ’| x') = 


2nh sin (f> 


x exp 


2Tl sin (f> . 


(x ' 2 + x 2 ) cos </> — 2xx' . (B40) 


6. Transformations generated by L z . Rotations 

Rotations are generated by A 12 in Sec. A 3. The trans¬ 
formation is given by 


U = exp 


* n 2 


(B41) 


The transformation rule for the energy operator is easily 
calculated by inserting a commutor and using (B3) as 
follows 


Upt& =Pt + U 


Pt 


U t 


= Pt + U\pt,t\ 


< 917 + 

~df 


= Pt + 9L Z . (B42) 


In order to obtain the transformation rules for the posi¬ 
tion and momentum operators we define 


X{0) = UxU\ 
Y{0) = UyU\ 


(B43) 

(B44) 


and calculate their derivatives with respect to the rota¬ 
tion angle 9 




x, L z 


y,L z 


t7+ =Y ( 0 ), (B45) 

C/t = —X (8 ). (B46) 


The solution to this system of differential equations to¬ 
gether with the boundary conditions X (0) = x and 
Y (0) = y yields the transformation rules 


UxW = cos Ox — sin By, 
UyU ] = sin&r + cos By. 


(B47) 

(B48) 


Following a similar procedure for the momentum opera¬ 
tors we obtain the rules 

UpJj X = cos 9p x — sin 9p y , (B49) 

Up v U ] = sin 9p x + cos Opy . (B50) 

The propagator can readily be obtained from the trans¬ 
formation rules (B47) and (B48) giving 


x,y 


17 + 


x', yj = S (x — x' cos 9 + y' sin 9) 

x 5 (y — x' sin# — y' cos#). (B51) 
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